Scalable experimental estimation of multipartite entanglement 
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Introduction.- Entanglement has been identified as a key 
resource for quantum information processing tasks. Further- 
more, it is clear that the dramatic advantage in using quan- 
tum mechanical systems instead of classical ones to process 
information emerges only in the limit of a very large num- 
ber of system components. Taking into account all the ex- 
tra qubits necessary for error correction, a quantum computer 
has to run on quantum registers composed of at least several 
thousand qubits to outperform its present-day classical coun- 
terpart. This explains the tremendous effort dedicated during 
the last few years to the experimental production and coher- 
ent manipulation of multiparticle entangled states of photons 
fl]|21|5]g||51|5), ions ||8ll9][T0l[TT1,andin cavity QED devices 
0. 

Also the experimental quantification of multipartite entan- 
glement has thus become a major issue of interest. In prin- 
ciple, such quantification can be carried out through quantum 
state tomography J^IEIEI, i.e., the complete reconstruction 
of the state's density matrix via the measurement of a com- 
plete set of observables, followed by the subsequent evalua- 
tion of a valid entanglement measure. In practice, however, to- 
mography rapidly saturates the available resources and is thus 
no viable strategy under the perspective of scalability. Clear 
evidence of this is given by the experimental characterization 
of genuine multiparticle entangled states of up to eight ions 
ifTTl : Ten hours of data aquisition - implementing measure- 
ments in 3 8 = 6561 detection bases, each corresponding to 
a different experimental setting - were followed by computa- 
tionally expensive data processing, to reconstruct the eight-ion 
density matrix of the experimentally prepared state. There- 
fore, full tomography of entangled ion chains composed of 
more than only eight ions appears largely impracticable. 

Quantum non-locality tests Q] 12 [4j [5] |6] [71 and entangle- 
ment witnesses B|8]|9l[T0][n][ll[T5][l6l provide alterna- 
tive means to assess the degree of entanglement of a quantum 
state, and were used in several experiments. Both these tech- 
niques require the measurement of only a few observables, but 
allow to detect the entanglement of only a small class of states. 
This implies that some a priori knowledge on the state to be 
analyzed is necessary. A simple entanglement measurement 
scheme for arbitrary mixed states is therefore highly desir- 
able. 

First steps in this direction were taken in IfTTl [T8ll . where 



multipartite concurrence [ 19 1 was shown to be directly acces- 
sible through projective measurements on two identically pre- 
pared quantum states. This original approach lfT71[T8ll - exper- 
imentally demonstrated for twin photon entanglement EOlETIl 
- was restricted to the ideal case of pure states, and a first gen- 
eralization for mixed states was given in ll22l . yet applicable 
only for bipartite systems. Here we come up with the ultimate 
formulation of this approach to direct experimental entangle- 
ment estimation, for mixed states of quantum systems with an 
arbitrary number of constituents. Our procedure, based on lo- 
cal parity measurements, features excellent scaling properties: 
the number of required observables - which can all be probed 
in one single experimental setting - is equal to the number of 
subsystems. 

The observable lower bound.- Consider two copies of an 
arbitrary mixed state g of an iV-partite quantum system with 
Hilbert space Ti. We introduce an observable V such that 

C 2 (g) >Tt[q®qV] , (1) 

i.e., that allows us to experimentally bound the concurrence 
of the state from below. The Hermitian operator V acts on the 
composite Hilbert space associated with the two-fold copy of 
the system fi^fi, and has the two following remarkable prop- 
erties: ( i) it can be detected through projective measurements 
of only N two-particle observables, and (ii) a single experi- 
mental setting is required throughout the detection process. 

The required two-particle measurements are simultaneous 
parity measurements on each particle and its copy. In each 
run of the experiment, measuring the local parity state of all 
N pairs defines an event in which each pair is projected onto 
either a symmetric or an antisymmetric state. From all possi- 
ble events we distinguish three types: 

(i) the entire system together with its copy is projected 
onto a globally symmetric state - symmetric with re- 
spect only to the exchange of both copies of the entire 
system; 

(ii) the entire system and copy are projected onto a globally 
antisymmetric state - antisymmetric with respect only 
to the exchange of both copies of the entire system; or, 

C Hi) system and copy are projected onto a full locally sym- 
metric state in which all N particle-copy pairs are si- 



2 




CD 



© 




(D 




© 


© 


BSM 


BSM 


BSM 


BSM 


BSM 


BSM 



© 



Figure 1: An eight-qubit state q, together with its copy, is en- 
coded, for example, in strings of two-level ions and is subject to 
Bell-state measurements (BSMs) on each pair. The entanglement 
of q is obtained from the joint probabilities of appearance of singlets 
and triplets. 



multaneously found to be in a symmetric state (which 
is a particular case of f ij). 

The probabilities of these three events suffice to obtain the 
expectation value of V, as described below. The measure- 
ment protocol is sketched in Fig. [T] for two strings of eight 
ions, reminiscent of the experimental situation in Hill . Each 
particle, together with its counterpart in the copy, is subject 
to a local parity measurement, which reduces to a Bell-state 
measurement, since the particles are qubits in this example. 
Recording the abundance of singlets in the string of eight ion 
pairs allows to infer the probabilities of the three above events 
immediately. 

From a more technical point of view, the system's Hilbert 
space H is a tensor product H = Hi Hi ® ... ® Ttjy of 
the single particle Hilbert spaces Hi, 1 < i < N. The sym- 
metric and antisymmetric subspaces Hi Hi and Hi A Hi 
of the Hilbert space Hi Hi of two copies of the i-th single- 
particle subsystem are defined as the subspaces spanned by all 
states that acquire a phase shift of or tt, respectively, upon 
exchange of the single-particle copies. These two subspaces 
are associated with the local two-particle projectors Pi and 
P l _ = 1 — P"\_. The globally symmetric and antisymmetric 
subspaces H H and H A H are the subspaces of all states 
that are symmetric and antisymmetric with respect to the ex- 
change of two copies of the entire system, and not only of 
some subsystems, and are in turn associated with the global 
projectors P + and P_ = 1 P + . In terms of these, our 
observable can be explicitly expressed as 



V = A P, 



(1 



,1-N 



)P- 



(2) 



Since the symmetric (antisymmetric) global projector P + 
(P_) can be decomposed into a sum of all products of N local 
projectors with an even (odd) number of antisymmetric local 
projectors, it suffices to measure the parity of the N pairs of 
copies to reconstruct Tr[g qV~\. 

Finally, it is important to note that V, as defined in has 
an equivalent interpretation to that of its bipartite analogue 
E2l . just with a much more intricate combinatorial struc- 



ture: the expectation value of A = 4 (P + 



pN 



g, however, a positive expectation value of A can have two 
causes: entanglement or mixedness of g. In turn, the oper- 
ator P quantifies the degree of mixing of g: 1 — Tr[g 2 ] = 
2Tr[g g P-]. The linear combination of A and P in (j2j) 
therefore rescales the expectation value of A with respect to 
the state's intrinsic impurity, and thus provides an estimate of 
the inscribed entanglement through a lower bound of multi- 
partite concurrence, as elaborated in the appendix. 

Tightness of the bound- Let us finally test the tightness of 
the observable bound on mixed random states. In Figs.[2]and[3] 
we plot the expectation value of the operator |2]), versus con- 
currence in quasi-pure approximation E3l (which is known 
to yield very good approximations for weakly mixed states), 
for 10 5 random states of 4-qutrit and 5-qubit systems, respec- 
tively, and for different degrees of mixing. Mixed states of dif- 
ferent purity were obtained by acting with the generalized de- 
polarizing channel (which essentially mixes a pure state with 
the identity) [24] onto 10 5 random pure states, for three dif- 
ferent coupling strengths. 

As spelled out by the comparison in Figs. [2] and [3] the ob- 
servable bound is hardly weaker than the quasi-pure approxi- 
mation. In fact, the comparison is excellent for weakly mixed 
or highly entangled states. On the other hand, for some very 
strongly mixed or very weakly entangled states other tech- 
niques involving few measurements, such as 'tailored wit- 
nesses' [25 1, may be used to improve the tightness of the en- 
tanglement estimation if some a priori knowledge of the state 
is available. The expectation value of however, provides 
a directly observable non-trivial bound for any unknown mul- 
tipartite state's concurrence. 

Conclusions .- We have derived a general lower bound for 
the entanglement of mixed quantum states, which provides a 
hierarchy of observable entanglement measures. As such, our 
result has the essential virtue of scalability for unknown, mul- 
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Figure 2: Observable lower bound versus concurrence, in quasi- 
pure approximation for 10 J four-qutrit density matrices with strong, 
intermediate, and weak mixing. The dashed line indicates equality 
of our present measurable bound and of entanglement in quasi-pure 
approximation. The tightness of the observable bound is excellent for 
strongly entangled or weakly mixed states, but remains surprisingly 
good even for strong mixing. 



yields the concurrence of pure states lfl8l . For a general state 
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Figure 3: Same as Fig. 2, now for 10 5 five-qubit density matrices at 
various levels of purity. Once more, the tightness of the observable 
bound is excellent for strongly entangled or weakly mixed states, and 
still surprisingly good for strongly mixed states. 



convex-roof decomposition, and inequality ([T} is automati- 
cally satisfied. Therefore, we seek V such that 

C^)C{</>) > <8 (4>\V\1>) 03 |</>) (3) 

holds for any two arbitrary pure states \ijj), \ <j>) G Ti. Such 
an observable is known for the bipartite concurrence c: v = 
4(P+ - Pj 03 P£ - |(P± 03 Pi + Pj 03 P^)) (221. Now, we 
can make use of the fact that the TV-partite concurrence can be 
decomposed into bipartite terms as 



cw^^/y^w, (4) 

where the sum is taken over the bipartite concurrencies Cj cor- 
responding to each subdivision of the entire system into two 
subsystems. This allows us to bound our quantity of interest 
from below as 



tipartite mixed quantum states in arbitrary finite dimensions. 
Given a two-fold copy of the state to be analyzed, our bounds 
are experimentally accessible, with linear scaling of the ex- 
perimental overhead with the number of system constituents. 
While derived for a specific type of multipartite concurrence 
|[T9l l26 l. equivalent expressions can be found for other observ- 
able multipartite concurrences with the same algebraic struc- 
ture ifTTl |26l . This defines a versatile toolbox for the exper- 
imental probing of quantum correlations inscribed into ever 
larger multicomponent quantum systems, an essential prereq- 
uisite for scaling up quantum information technology. 

Appendix- Here we prove that the observable defined in (|2| 
satisfies Eq. ([TJ, for any state g: The concurrence of g is given 
by the convex roof ETl C(g) = inf J^j Cfyj)' i- e - me min- 
imal average concurrence over all (subnormalized) pure-state 
decompositions g = £\ \^j)(tpj\- If Ejfe C '(V , i)C('0fc) > 
Tr[g <g> q V] = £jfc<V»j| (3 (^fe|V|Vj> <8> |V»fc> holds for 
all decompositions {\ipj)}, then it also holds for the optimal 



C(*)C?(*) = 2 2 - w ^ c fW^ c 2(ci>) (5) 

> 2 2 - N J2^)cm (6) 

i 

> 2 2 - JV ^<V|03(^| Vi |V)®k) , (7) 



w here w e made use of the Cauchy-Schwarz inequality 

\/Hi x \ VJ2iVi > J2i x iVi and the above knowledge on 
bipartite systems. It is now a matter of straightforward alge- 
braic gymnastics to show that V — J2i v ^ wna t finishes the 
proof of Eq. 
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